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Abstract 

Heterodyne holography is an extremely versatile and powerful 
holographic technique that is able to fully control the amplitude, 
phase and frequency of both the illumination and holographic ref¬ 
erence beams. Applied to vibration analysis, this technique is 
able to detect the signal of the carrier or the one of any vibration 
sideband. A complete analysis of the vibration signal can thus 
be made, and 2D map of the vibration amplitude and phase in 
all points of the surface of the vibrating object can be obtained. 

Since the sensitivity is limited by shot noise, extremely low vibra¬ 
tion amplitude (< 0.01 nm) can be studied. 
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1 Introduction 

Digital holography is a fast-growing research field that has drawn 
increasing attention. The main advantage of digital holography is 
that, contrary to holography with photographic plates, the holo¬ 
grams are recorded by a CCD, and the image of the object is 
digitally reconstructed by a computer, avoiding photographic pro¬ 
cessing pp. To extract two quadratures of the holographic signal 
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(i.e. to get the amplitude and the phase of the optical signal) 
two main optical configurations have been developed: off-axis and 
phase-shifting. 

Off-axis holography is the oldest and the simplest configuration 
adapted to digital holography. In that configuration, the reference 
or local oscillator (LO) beam is angularly tilted with respect to 
the object observation axis. It is then possible to record, with a 
single hologram, the two quadratures of the object complex field 
[T]. However, the object field of view is reduced, since one must 
avoid the overlapping of the image with the conjugate image alias. 
Phase-shifting digital holography [2] makes possible to get phase 
information on the whole camera area by recording several images 
with a different phase for the reference (called here local oscillator 
or LO) beam. It is then possible to obtain the two quadratures 
of the field in an on-axis configuration even though the conjugate 
image alias and the true image overlap, because aliases can be 
removed by taking image differences. In a typical phase-shifting 
holographic setup, the phase of the reference is shifted by moving 
a mirror with a PZT. 

On the other hand, there is a big demand for full field vibra¬ 
tion measurements, in particular in industry. Different holographic 
techniques are able to image and analyze such vibrations. Double 
pulse holography 011] records a double-exposure hologram with 
time separation in the 1...1000 /iS range, and measures the instan¬ 
taneous velocity of a vibrating object from the phase difference. 
The method requires a quite costly double pulse ruby laser system, 
whose repetition rate is low. Multi pulse holography 0 is able to 
analyse transient motions, but the setup is still heavier (4 pulses 
laser, three cameras). 

The development of fast CMOS camera makes possible to an¬ 
alyze vibration efficiently by triggering the camera on the motion 
in order to record a sequence of holograms that allows to track the 
vibration of the object as a function of the time [Hill- The analysis 
of the motion can be done by phase difference or by Fourier anal¬ 
ysis in the time domain. The method requires a CMOS camera, 
which can be costly. It is also limited to low frequency vibra¬ 
tions, since a complete analysis of the motion requires a camera 
frame rate higher than the vibration frequency, because the band¬ 
width BW of the holographic signal, which is sampled at the cam¬ 
era at angular frequency ojccd must be lower than corresponding 
Nyquist-Shannon limit: BW < ^ ojccd- 

For a periodic vibration motion, the bandwidth BW is close 
to zero. Measurements can thus be done with much slower cam¬ 
eras. Powell and Stetson [S] have shown for example that an har¬ 
monically vibrating object yields alternate dark and bright fringes, 
whose analysis yields informations on the vibration motion. Picard 
et al. 0 has simplified the processing of the data by performing 
time averaged holography with a digital CCD camera. Time av¬ 
eraged holography has no limit in vibration frequency and do not 
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involve costly laser system, nor an expensive CMOS fast camera. 
Although the time-averaging method gives a way to determine the 
amplitude of vibration m quantitative measurement remain quite 
difficult for low and high vibration amplitudes. 

Heterodyne holography [mna is a variant of phase shifting 
holography, in which the frequency, phase and amplitude of both 
reference and signal signal beam are controlled by acousto optic 
modulators (AOM). Heterodyne holography is thus extremely ver¬ 
satile. By shifting the frequency ujlo of the local oscillator beam 
with respect to the frequency ujq of illumination, it is for example 
possible to detect the holographic signal at a frequency uj different 
than illumination loq. This ability will be extremely useful to ana¬ 
lyze vibration, since heterodyne holography can detect selectively 
the signal that is scattered by the object on a vibration sideband 
of frequency uJm = wq -I- mujA, where uja is the vibration frequency 
and m and integer index. 

In this chapter we will first present in section [5] heterodyne 
holography and its advantages in section 12.11 and 12.21 Then in 
section |31 we will apply heterodyne holography to vibration analy¬ 
sis. We will show in section how heterodyne holography can be 
used to detect the vibration sidebands, in section [23] how this side¬ 
band holography can be combined with stroboscopic techniques to 
record instantaneous velocity maps of motion, and in sections 13.dl 
and 13.51 how it can retrieve both small and large vibration ampli¬ 
tudes. 


2 Heterodyne holography 

Let us first describe heterodyne holography. A example of setup is 
shown on Fig. (Tj The object is an U.S. Air Force (USAF) resolu¬ 
tion target whose hologram is recorded in transmission geometry 
(the target is back illuminated). The camera records the interfer¬ 
ence of the signal field £{t) of optical angular frequency wq, with 
the reference (or local oscillator field) £Lo{t) of optical frequency 

wlo = Wo - Aw (1) 

where Aw is the frequency shift. Let us introduce the slowly vary¬ 
ing complex fields E and E^o- 

£{t) = -k (2) 

£Lo{t) = ELoe^‘^^°* + 

where E* and E^q are the complex conjugates of E and and 
j is the imaginary unit. The camera signal I is proportional to the 
intensity intensity of the total field \£{t)+£Lo{t)\^- We have thus: 

I = jAio -k A (3) 

= IAloI" + \E\^ + EloE + ElqE* 
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Figure 1: Typical heterodyne holography setup. L: main laser; BS: beam 
splitter; AOMl, AOM2: acousto-optic modulators; BE: beam expander; 
M, mirror; Al, A2: light attenuators; USAF, transmission USAF 1951 
test pattern; CCD: camera; El, Ej, Elo: E: laser, illumination, 

reference (i.e. local oscillator LO) and object fields; uiaomi/ 2 - driving 
frequencies of the acousto optics modulators AOMl and AOM2; 6: off- 
axis angular tilt. 


Since the camera signal is slowly varying, we have neglected in Eq. 
[2]the fast varying terms (which evolve at frequency ~ 2a;o). More¬ 
over, to simplify the present discussion, we have not considered 
the spatial variations of I, E and Elo with x and y, in particular 
the spatial variations that are related to the off axis tilt angle 0 
of Fig. [T] In equation O E^qE is the -1-1 grating order 

term, that contains the useful information (since this term is pro¬ 
portional to E). The others terms: {EloI^ + (zero grating 
order term), and ElqE* ( -1 grating order term or twin 

image term) are unwanted terms that must be cancelled. To filter 
off these terms, and to select the wanted -1-1 grating order signal, 
4-phase detection is made. The AOMs driving angular frequencies 
waomi and ujaom 2 are tuned to have: 

Auj = UJAOM2 — ^^AOMl = I^CCd/^ (4) 

where loccd is is the angular frequency of the camera frame rate. 
Lets us consider 4 successive camera frames: Iq, A ••■A that are 
recorded at times t = 0, T ... 3T with T = 2Tr /{ujccd)- For these 
frames, the phase factor is equal to 1, j, -1 and —j. We 
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thus get: 


h = \Elo?+\E\^ + EI 0 E+ EloE* (5) 

h = \Elo?+\E\^+ j EloE-j EloE* 
h = \Elo?+ \E\^ - EloE- EloE* 
h = \ELof+\E\^-j EloE + j EloE* 

By linear combination of 4 frames, we get the 4-phase hologram 
E[ that obeys the demodulation equation: 

H = {lo - I 2 ) + j{h - h) (6) 

= ^EloE 

As wanted, the 4-phase demodulation equation yields a quantity 
H that proportional to the signal complex field E. Note that the 
coefficient AElq is supposed to be known and do not depends on 
the object. 

Heterodyne holography exhibits several advantages with re¬ 
spect to other holographic techniques: 

1. The phase shift is very accurate; 

2. The holographic detection is shot noise limited; 

3. Since the holographic detection is made somewhere near wlo 
(depending on the demodulation equation that is chosen), it 
is possible to perform the holographic detection with any fre¬ 
quency shift with respect to the object illumination angular 
frequency ujq. 

The first advantage will be discussed in section 12.11 the second 
one in section while the third one, which is the heart of the 
sideband holographic technique used to analyse vibration, will be 
discussed in section o 

2.1 Accurate phase shift and holographic detec¬ 
tion bandwidth 

In the typical heterodyne holography setup of FiglU the signals 
that drive the AOMs are generated by frequency synthesizers, 
phase-locked with a common 10 MHz clock.. The phase shift 
Atf = AojT (that is equal to tt/2 in the 4 phase demodulation 
case) can thus be adjusted with quartz accuracy. 

To illustrate this accurate phase shift, holograms of a USAF 
target have been recorded with the Fig [T| setup, while sweeping the 
LO frequency ujlo- Figure [2] shows USAF reconstructed images 
that are obtained for different values of the frequency shift loq—lolo 

m- 

1 . For ujq — ujlo = wcc£)/ 4 =4-2.5 Hz, the image that is recon¬ 
structed is sharp and correspond to the 4-1 grating order. On 
the other hand, the -1 grating order signal is very low; the 
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Figure 2: USAF target reconstructed images (a,b) displayed with linear 
grey scale for the reconstructed field intensity \E\^ with ojq — ojlo = 
wccd/ 4 = +2.5 Hz (a) and wq - uolo = -ujccd/^ = -2.5 Hz (b). 
Zooms (c,d) of the bright zones of the (a) and (b) images. Images (a,b) 
and zoom (c,d) correspond to the m = +1 (a,c) and m = —1 (b,d) 
grating orders. 


magnitude of the -1 parasitic twin image image is negligible 
in front of the +1 contribution. Image and zoom of image 
are seen on Figl2](a) and (c). 

2. For Wo — wlo = —wcco/4 = —2.5 Hz, the image that is 
reconstructed is blurred and corresponds to the -1 grating 
order twin image. The magnitude of the +1 grating order 
image is negligible in front of its -1 counterpart. Image and 
zoom of image are seen on Figl^(b) and (d). 

We have measured the total field energy 17± in both +1 and -1 
grating order reconstructed images: 

U± = (7) 

x,V 

where y pixels of the reconstructed field 

image E{x,y) corresponding to the ±1 grating order zones. As 
seen on Fig. [21 the energy of the signal that is measured in the 
+ 1 grating order (i.e. 17+) is maximum for ujq — lolo = ^cgd!^ = 
+2.5 Hz, and null wq — w+o = —2.5 Hz. Similarly, C7_ is maximum 
for Wo — w+o = —2.5 Hz, and null for wq — w+o = +2.5 Hz. By 
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Figure 3: Field energy of the +1 (a) and -1 (b) grating order. The 
vertical axis is in logarithmic scale arbitrary units. The points corre¬ 
spond to the experimental data obtained by sweeping the LO frequency 
with (wq — ojlo) = to -|-5 Hz with 0.1 Hz increments. The solid gray 
curves are +U+{uj — ojlo) (a) and +U-{uj — ojlo) (b) given by Eql5]in 
the same frequency range: (cu — ojlo) = to -|-5 Hz. 


adjusting wq — ojlOi h is then possible to select the grating order 
that is detected. 

For 4 phase detection with a local oscillator of frequency ujlo, 
the total energy U± = \E\^ detected at frequency uj in the ±1 
grating order can be easily calculated M- 


U±{lj — ujlo) 


n—3 


AT' 


E(±j)' 


nnT+T' 


dt 


n=0 


' t—nT 


( 8 ) 


sinc(7r(a; - uilo)T') ■n^jn{u,-u,Lo)T 
4 


E/ 

n=0 


where sinc(a;) = sinx/x. In equation [S] T is the frame period, 
and T' the exposure time. The coefficient 1/4T' is thus a nor¬ 
malization factor. On the other hand, the factors (±j)" corre¬ 
sponds to the coefficients of the demodulation equation (see Eq. 
O, and jg interference instantaneous phase factor, 

which must be integrated over the exposure time from t = nT to 
t = nT + T'. 

In the experiment of FigO we have T = T' = 2-k/ujccd = 100 
ms. We have plotted the field energy U± given by Eql5]as a func¬ 
tion of a; — ujlo on Fig. [3] (solid grey lines). As seen, experi¬ 
ment agrees with the theoretical curve of Eq. [Sj The shape of 
the curves represents here the frequency response spectrum of the 
holographic device considered as a detector. For the -1-1 grating 
order, detection is centered at frequency w = ljlo + ojccd/^- The 
measurement bandwidth BW is 2.5 Hz. It is equal to the inverse of 
the measurement time of 4 frames i.e. BW = l/4r. It illustrates 
the coherent character (in time) of holographic detection. 
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Figure 4: Theoretical field energy U+{uj — ojlo) given by EqE] plotted 
in linear (a) and logarithmic scale (b) as a function of {u — ojlo) in Hz 
Units for the +1 grating order. The number of frame is Umax = 4 (solid 
grey line), Umax = 8 (solid black line) and Umax = 16 (dashed black 
line). 


It is possible to increase the selectivity of the holographic co¬ 
herent detection by increasing the measurement time, i.e. by in¬ 
creasing the number of frames Umax used for demodulation. In 
that case, Eq. and EqlS] must be replaced by similar equations 
involving Umax frames in place of 4. We have: 

rijnax — ^ 

£ fin (9) 

n=0 


and 

U±{uj — iOLo) 


1 


rimax — ^ 


pnT -\-T' 

(±jT / dt 


■II 

n=0 

sinc(7r(a; — u}lo)T') 


J t—nT 

n-max — ^ 

E 

n—O 


n^jn{uj-uJLo)T 


Figure m plots the detection frequency spectrum U+ for Umax = 4, 
8 and 16. As seen, the detection bandwidth BW decreases with 
nmax- It is equal to: 


BW = 


1 


'^max 


T 


(11) 


The noise, uniformly distributed in frequency, is expected to de¬ 
crease accordingly. We will show in section 12.21 that very high de¬ 
tection sensitivity can be obtained by holographic detection with 
large number of frames Umax- 


2.2 Shot noise holographic detection 

Because the holographic signal results in the interference of the 
object signal complex field E with a reference (LO) complex field 
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Figure 5: Heterodyne detection of 1 photon per pixel of signal with N 
photon local oscillator. 


Elo whose amplitude can be much larger (i.e., E^o ^ E), the 
holographic detection benefits from ’’heterodyne” or ’’holographic” 
gain (i.e., IEE'Iq] and is thus well suited for detection 

of weak signal fields E. Holographic detection can reach the the¬ 
oretical limit of noise which corresponds to a noise equivalent sig¬ 
nal of 1 photo electron per pixel during the total measurement 

time [isiiiniiiiiizi. 

To illustrate this point, let us consider the interference of a 
weak signal of about 1 photo-electron per pixel and per camera 
frame, with a large local oscillator reference signal of IV = 10^ 
photo-electrons: see Fig. |5] In this example, the signal intensity 
is \E\'^ = 1, while the holographic interference term {EE’IqI = 100 
is much larger. The detected signal is \Elo + E\'^ ~ 10^ photo¬ 
electrons. 

Because of the quantum nature of the process involved in dig¬ 
ital holography (laser emission, photodetection...), the detected 
signal in photo electron Units is random Gaussian integer number, 
whose average value is N = 10“^, and whose standard deviation is 
'/N = 100. These fluctuations of the number of photo electrons 
are called shot noise. Here, with 1 photo electron of signal, the 
heterodyne signal is 100 photo electrons and the noise 100 too. 
The noise equivalent signal (for the energy |ifp) is thus 1 photo 
electron per pixel and per frame. 

Let us now study how shot noise varies with the number of 
frames Umax used for detection. As in any detection process, the 
noise in energy is proportional to the measurement time and to the 
detection bandwidth BW. Since the measurement time is UmaxT, 
and since the detection bandwidth is BW =1 /rimaxT (see Eq. [TT] 
), the shot noise in energy do not depends on rimax ■ The shot noise 
equivalent signal is thus the same than for one frames. It remains 
equal to 1 photo electron per pixel whatever the number of frames 
"r^max and measurement time UmaxT are. 

Let us now discuss the ability to reach this shot noise optimal 
sensitivity in real life holographic experiments. Since we consider 
implicitly a very weak signal, the noises that must be considered 
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Figure 6: Fourier space hologram intensity (i.e. \H{kx-,kyf‘) displayed 
in arbitrary log scale for one phase {H = Iq) and 4 phases {H = (Iq — 
h) + Jih — I3) detection. 


are 

1. the read noise and dark current of the camera, 

2. the quantization noise of the camera A/D converter, 

3. the technical noise of the LO beam, 

4. and the LO beam shot noise, which yield the theoretical noise 
limit. 

For a typical camera, the full well capacity is 2 10"^ photo electrons, 
and a good practice is to work with TV = 10"^ photo electrons for 
the local oscillator. Shot noise on the camera signal (100 photo 
electrons) is thus much larger than the camera read noise (1 to 20 
photo electrons) and the camera dark current (a few photo elec¬ 
trons per second). If the camera is 12 bit, the full well capacity 
corresponds ^ 2^^ digital count (DC). As a results, the quantiza¬ 
tion noise (^ 7 photo electrons) can be neglected too. 

Let us discuss on the technical noise of the LO beam, and on 
the way to filter it off. 

1. We have displayed, on Fig. 0 examples of Fourier space holo¬ 
gram H{kx,ky) = FFT( H{x,y)) (where FFT is the Fourier 
transform operator). The signal fields yields the interfer¬ 
ence term EE^q that is located in the left hand side of the 
Fourier space images of Fig. El (a) and (b). It corresponds 
to the grating order -|-1. The local oscillator signal \Elo\^ 
is located in the center of the Fourier space. It is visible on 
Fig. El (a) that is obtained with 1 phase hologram. As seen, 
the LO signal are separated in the Fourier space because of 
the off axis configuration. It is then possible to filter off the 
parasitic LO signal by a proper spatial filtering in the Fourier 
space as shown by m- As noticed, this operation filters-off 
both the LO signal and the LO technical noise. 

2. On the other hand, because of phase shifting, the signal and 
LO fields E and E^o evolve at two different time frequencies 
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Figure 7: Reconstructed image of a USAF target in dim light. Hologram 
are recorded by 4 phases heterodyne holography with T = 100 ms, and 
nmax = 600 (a,b) and Umax = 6000 (c). The coherent measurement 
time UmaxT is thus 1 minute (a,b) and 10 minutes (c). Illumination is 
adjusted so that the USAF signal integrated is 1 photon per pixel over 
the whole measurement time UmaxT in (a). Holograms (b) and (c) are 
recorded by adding a neutral density filter D=1.0 on illumination. The 
USAF signal is thus 0.1 photon per pixel in (b) and 1 photon per pixel 
in (c). Display is made with arbitrary linear scale for intensity. 


u) and ujlo with oj = ujlo + ^ccd/^- It is then possible to 
filter off the LO signal in time. This is done by the 4 phases 
demodulation process since the LO term IAloP vanish in 
Eq. and [Till This filtering in time is illustrated by Fig. El 
since the 0 grating order signal that is large for single phase 
off-axis holography (Fig. [S](a)), roughly vanishes in the four 
phase case (Fig. E](b)). Time filtering is also illustrated by 
FigOand Fig|3]and Eq. Eland [101 since the detected energy 
U+/_ vanishes for u — ulo = 0 i.e. for detection at local 
oscillator frequency ujlo- Here again, one filters off both the 
LO signal and the LO technical noise. 

By combining off axis and phase shifting, it is then possible apply 
to a double filtering (in space and time) that filter off the LO 
technical noise very effectively. One can gets then the shot noise 
ultimate sensitivity that corresponds to a noise equivalent signal 
of 1 photon per pixel whatever the measurement times is. 

Figure|3illustrates heterodyne holography ability of imaging an 
object (here an USAF target) in dim light with shot noise limited 
sensitivity. In figure [71 (a), the measurement time is 1 minute and 
the illumination power is adjusted such a way the USAF signal is 
about 1 photon per pixel in 1 minute. The visual quality of USAF 
image is quite good. We can say that SNR is about 1. In figure |7| 
(b), the measurement time is the same, but illumination is divided 
by 10 by using a neutral filter of density D = 1.0. SNR is low, 
and one cannot see the grooves of the USAF target. Figure [7] (c) 
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Figure 8: Heterodyne holography setup applied to analyse vibration of a 
clarinet reed. L: main laser; AOMl, AOM2: acousto-optic modulators; 
M: mirror; BS: beam splitter; BE: beam expander; CCD: camera; LS: 
loud-speaker exiting the vibrating clarinet reed at frequency 


is obtained with the same illumination level that Fig. [3 (b) (i.e. 
with neutral density filter) but the measurement time is multiplied 
by 10 (i.e. 10 minutes in place of 1 minute). The visual quality 
of reconstructed image is the same than for Fig. [3 (a) with SNR 
about 1. This experiment shows that the reconstructed image 
quality depends on the total amount of signal, and do not depends 
on the time needed to get that amount of signal: to get SNR ^ 1 
one needs about 1 photon per pixel, for any measurement time. 


3 Holographic vibrometry 

Lets us now apply heterodyne holography to vibration analysis. 

3.1 Optical signal scattered by a vibrating ob¬ 
ject 

Consider a point of the objet (for example a clarinet reed) that 
is studied by heterodyne holography (see Fig. [5]), vibrating at 
frequency oja with amplitude Zmax- The displacement z along the 
out of plane direction is 

z{t) = ZmaxSinUJAt (12) 

In backscattering geometry, this corresponds to a phase modula¬ 
tion (fit) of the signal: 

ip{t) = 47rz(t)/A (13) 

= <i)sina;,4t 
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Sideband index rn Sideband index rn 

Figure 9: Relative amplitude \Em\‘^/\E\2 = of the sideband 

component m for an amplitude modulation of the phase equal to ‘h = 0.3 
(a) and = 3.0 (b) rad. 


where A is the optical wavelength and d* the amplitude of the phase 
modulation of the signal at angular frequency oja- 

= 47rZmaa:/A (14) 

Let us define the slowly varying complex amplitude E{t) of the 
field £{t) scattered by the vibrating object. We have: 

£{t) = + c.c. (15) 

Because of the Jacobi-Anger expansion, we get: 

E{t) = = A 

= R ^ J™(4>) 

m 

where E is the complex amplitude without vibration, and Jm 
the mth-order Bessel function of the first kind, with J-m{z) = 
— for integer m and real 2 . The scattered field £{t) is 

then the sum of the carrier and sideband field components £m{t) 
of frequency Wm, where m is the sideband index with: 

+ 00 

£{t) = Y. ^rn{t) (16) 

m— — oo 

£m{t) = 

LOm = Wo + mu: A 

where E^ is the complex amplitude of the field component £m(t). 
Equation [T^ yields: 

Em = Jmm E (17) 

Figure [2] presents the distribution of the field energy on the 
sidebands components If the amplitude of modulation d* 
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Figure 10: Reconstructed holographic images of a clarinet reed vibrating 
at frequency = 2143 Hz perpendicularly to the plane of the figure. 

Figure (a) shows the carrier image obtained for m = 0. Fig. (b)-(d) 
show the frequency sideband images respectively for m = 1, m = 10, 
and m = 100. A logarithmic grey scale has been used 


is low, most of the energy is on the carrier |ifo|/|Ap — !> and 
energy |ifmP decreases rapidly with the sideband index m. If the 
amplitude 4) is large, the energy of the carrier is low lAol/lAp 1, 
while energy is distributed over many sidebands ■ 

3.2 Selective detection of the sideband compo¬ 
nents Em- sideband holography 

Heterodyne holography that is able to perform the holographic de¬ 
tection with any frequency shift uj — ujq with respect to the object 
illumination angular frequency wq is well suited to detect the vi¬ 
bration sideband components Em ■ To selectively detect by 4 phase 
demodulation the sideband m of frequency ojm, the frequency ujlo 
must be adjusted to fulfil the condition : 

ULO = ^m—^CCDl^ (18) 

= Wo + ruwA - wccn/4 

Figure [12] shows images obtained by detecting different side¬ 
band TO of a clarinet reed m- The clarinet reed is attached to a 
clarinet mouthpiece and its vibration is driven by a sound wave 
propagating inside the mouthpiece, as in playing conditions, but 
the sound wave is created by a loudspeaker excited at frequency 
and has a lower intensity than inside a clarinet. The excitation 
frequency is adjusted to be resonant with the first flexion mode 
(2143 Hz) of the reed. 

Figure [Tni(a) is obtained at the unshifted carrier frequency ujq. 
It corresponds to an image obtained by time averaging holography 
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[S]. The left side of the reed is attached to the mouthpiece, and the 
amplitude of vibration is larger at the tip of the reed on the right 
side; in this region the fringes become closer and closer and difficult 
to count. The mouthpiece is also visible, but without fringes since 
it does not vibrate. Similar images of clarinet reeds have been 
obtained in [201 [H], with more conventional techniques and lower 
image quality. Figures (TO] (b) to (d) show images obtained for the 
sidebands m = 1, 10 and 100. As expected, the non-vibrating 
mouthpiece is no longer visible. Figure [TO] (b) shows the m = 
1 sideband image, with Ji fringes that are slightly shifted with 
respect to those of Jq. Figure [TO] (c) shows the image of sideband 
m = 10 and Fig. [TO] (d) the m = 100 one. The left side region 
of the image remains dark because, in that region, the vibration 
amplitude is not sufficient to generate these sidebands, Jm{z) being 
evanescent for z < m. 

To quantitatively visualize the vibration amplitude cuts of 
the reconstructed images signal \Em{x,y)\’^ are made for differ¬ 
ent sideband orders m along the horizontal line y = 174. This 
value has been chosen because it corresponds to a region where 
the fringes are orthogonal to the y axis. To build the central part 
of Fig. [TT] the position of the antinodes of \Em{x, 174)^ are put in 
correspondence with the antinodes of |Jm(4>)p. Correspondence 
is made for m = 0, 1, 5, 10, . . . 100. Note that this method is 
insensitive to inhomogeneities in the illumination zone. Therefore, 
no normalization is required. The curve seen in the central zone of 
Fig. [TT] represents the amplitude of phase oscillation 4) in radian 
as a function of pixel index along x for y = 174. 

Remark: In a typical heterodyne holography setup, digital 
synthesizers drive the acousto optic modulator at ujaomi and 
^AOM2, the camera at locCD, and the vibration frequency at loa- 
These synthesizers use a common 10 MHz reference frequency, and 
are driven by the computer. It is then possible to automatically 
sweep ujAi and loaomi (or WAOM2) in order to fulfil Eq. [THl so 
that detection remains ever tuned on a given sideband. Figure [T2] 
shows an example m- A series of 26 x 7 images of a clarinet reed 
are obtained on sideband m = 1 by sweeping the frequency loa 
from 1.4 kHz up to 20 kHz by steps of 25 cents. The amplitude 
of the excitation signal is exponentially increased in the range 1.4 
to 4 kHz, from 0.5 to 16 V, then kept constant at 16 V up to 20 
kHz. This crescendo limits the amplitude of vibration of the first 
two resonances of the reed. The different vibration modes of the 
reed can be easily recognized on the reconstructed reed images of 
Fig. [T2] 

3.3 Sideband holography for large amplitude of 
vibration 

In the previous section (section 15.21 and Fig. fTTT) we have shown 
how the comparison of dark fringes for different sideband leads to 
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Figure 11: A slice of the data along the y = 174 line is used in this 
figure; the x horizontal axis gives the pixel index (100 pixels correspond 
to 3.68 mm.), the vertical axis the vibration amplitude 4>. The lower 
part of the figure shows the normalized signal |£'m(x)p/|£'(x)p (where 
E(x) is obtained loud-speaker off) for m = 5, with a downwards axis. 
The left part shows the corresponding square of the Bessel function 
|J5(4>)P with a leftward axis. The zeroes of the two curves are put in 
correspondence, which provides the points in the central figure. Similar 
correspondences are made harmonic order m = 0,1,5,10...100. Different 
gray densities are used for different m. The crosses correspond to m = 5. 
The juxtaposition of the points for all values of m gives an accurate 
representation of the amplitude of vibration as a function of x. 
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Figure 12: Clarinet reed reconstructed images obtained on sideband 
m = 1. Frequency uja is swept from 1.4 KHz up to 20 KHz, and images 
are displayed from left to right and top to bottom (26 x 7 images). The 
display is made with arbitrary grey scale for the intensity |Ko(x,y)p. 


a determination of the vibration amplitude 4>(a;,y) at each point 
of the object. This determination is non-local, since it involves 
counting fringes from one reference point of the image to the point 
of interest, so that large amplitudes are not accessible. The vi¬ 
bration amplitude 4)(a:, y) can be determined by another approach 
that completely eliminates the necessity of counting fringes, giv¬ 
ing a local measurement of the amplitude of vibration <i>, even for 
large 4) [23] . 

For large amplitude of vibration ( 4) ^ 1), the distribution 
of the sideband energy over m exhibits a sharp variation 

from maximum to zero near m ~ 4>, as seen on Fig. [T3] that 
plot \Em\^/\E\^ for 4) = 30.3. This property can be understood 
if one consider the limits Zmax ^ A. In that case, one can define 
an instantaneous Doppler angular frequency shift ujoit), and an 
instantaneous sideband index m{t) that are continuous variables : 

WD(t) = a;A4> cos(a;At) (19) 

m{t) = = (1) cos(a;At) 

WA 

Because of its sinusoidal variation, m{t) spend more time near the 
extreme points m = ±4>. The Doppler continuous distribution 
spectrum of m{t) that is displayed in light grey shade on Fig. [13] 
are thus maximum near the extreme points m = ±4>. 

The vibration amplitude 4>(a:,?/) can be determined for each 
location x,y^ by measuring \Em{x,y)\'^ for all sidebands m, and 
by determining for each location cc, y the sideband index m that 
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Sideband index m 


Figure 13: Relative intensity /\E\2 = of the sideband 

component m for an amplitude modulation of the phase equal to = 
30.3 radiant. The light grey shade shows the Doppler spectrum obtained 
from the vibration velocity distribution, with a continuous variable m. 
Both spectra fall abruptly beyond m = 30.3, which corresponds to the 
Doppler shift associated with the maximum velocity. 
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Figure 14: Cube of data obtained from the reconstructed holographic 
images of a vibrating clarinet reed. The sideband images with m = 0, 
20, 40 ...120 are shown in arbitrary linear scale. The white dashed lines 
correspond to x = 249 and y = 750, i.e. to the point chosen for Fig. [T5J 


correspond to a fast drop of the signal \Ejn{x,y)\^. The method 
is robust and can easily be used even when the fringes become so 
narrow that they cannot be resolved, which gives immediate access 
to large amplitudes of vibration. 

By successively adjusting the frequency ojlo of the local os¬ 
cillator to appropriate values, one records the intensity images 
\Em{x, y)\^ of the sidebands as a function oi x y and m. One then 
obtains a cube of data with three axes x, y and m. Figure [T?] shows 
the images obtained for for to = 0, 20 ...120 that correspond to 
cuts of the cube along x, y planes. The images illustrate how, when 
TO increases, the fringes move towards regions with larger ampli¬ 
tudes of vibrations. Since the right part of the reed {x > 800) is 
clamped on the mouthpiece, no signal is obtained in regions near 
the clamp where 4> = Anzmax/^ < xn (right part of the images). 

Figurefini'a'l displays a 2D cut (coordinates x, to) of the cube of 
data along the horizontal plane y = 750 (horizontal white dashed 
line in Fig. [T^ . The envelope of the non-zero (non black) part of 
the image provides a measurement of the amplitude of vibration 
in units of A/47r. One actually obtains a direct visualization of 
the shape of the reed at maximal elongation, from the right part 
clamped on the mouthpiece to the tip on the left. The maximum 
amplitude correspond to 4) = 120 rad. 

Figure [T^ shows images obtained at higher excitation ampli- 
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Figure 15: Images corresponding to cuts of 3D data of the reconstructed 
images along the planes y = 750 (a) and x = 249 (b). Fig (a) shows 
the deformation of the object along its axis and Fig (b) a transverse cut 
with a slight vibration asymmetry. A logarithmic intensity scale is used. 


512 

y 

1023 
1250 

n 

-1250 

Figure 16: (a): image reconstructed with sideband n = 330, with a large 
amplitude of vibration, (b) is the equivalent of Fig. fT^ai. but with 
positive and negative m values. One measures a maximum vibration 
amplitude of Zmax — 60 ym. A logarithmic intensity scale is used. 
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tudes, about 10 times larger that for FiglT^ Figure [T^ (a) shows 
the images obtained for m = 330: the fringes are now com¬ 
pletely unresolved, but the transition from zero to non-zero in¬ 
tensity remains very clear. With a single hologram, and with¬ 
out fringe counting, one obtains a clear marker of the line where 
y) = 330 rad. Figure fW bl shows the equivalent of Fig. fTM ab 
but with a higher excitation level, and this time for positive and 
negative values of m. Data range up to about \m\ ~ 1140, cor¬ 
responding to Zmax — 58.4 yum. Since the vibration amplitude is 
much larger than A, the continuous approximation for m is valid, 
and the images of Fig. [T^can be reinterpreted in term of classical 
Doppler effet. 

Taking advantage of the sideband order m of the light scattered 
by a vibrating object adds a new dimension to digital holography. 
Each pixel x^y ol the image then provides an information that is 
completely independent from others, which results in redundancy 
and robustness of the measurements. Looking at the edges of the 
spectrum provides an accurate determination of the vibration am¬ 
plitude and avoids a cumbersome analysis of the whole cube of data 
cube, giving easy access to a measurements of large amplitudes of 
vibration. 

3.4 Sideband holography with strobe illumina¬ 
tion 

Both time averaged [U] and sideband digital holography m record 
the holographic signal over a large number a vibration periods. 
These two techniques are not sensitive to the phase of the vibra¬ 
tion, and are thus unable to measure the instantaneous velocities 
of the object. To respond this problem, Leval et al. m com¬ 
bine time averaged holography with stroboscopic illumination, but, 
since Leval uses a mechanical stroboscopic device, the Leval tech¬ 
nique suffer of a quite low duty cycle (1/144), and is limited in low 
vibration frequencies (a;^/27r < 5 kHz). 

To overcome these two limitations, it is possible to combine 
sideband digital holography with stroboscopic illumination and de¬ 
tection synchronized with the vibration [25] . This can be achieved 
without changing the experimental sideband holography setup of 
Fig [5] by switching electronically on and off the Radio Frequency 
signals that drive the AOMs at ujaomi and ujaom2 — 80 MHz. 
Figure [T7] shows a typical chronogram of stroboscopic illumination 
and detection. 

Without stroboscopic illumination (and detection), the Doppler 
velocity spectrum is similar to the ones observed in Fig [12] section 
o It covers the entire range of speeds that can be reached during 
the sinusoidal motion, that is 4> < to < +4i, with two maxima near 
TO ~ ±4>. With stroboscopic illumination, the Doppler spectrum is 
modified, since it corresponds to the narrower velocities range that 
is reached during gated illumination. On the other hand, the finite 
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Figure 17: Chronogram of the signals. Sinusoidal signal of period Ta = 
2'kIuja exciting the reed (a). Rectangular gate that is applied to both 
illumination and reference beams (b) .Gate delay is with respect to the 
origin of reed sinusoidal motion. Gate duration is bt. Gated AOMl and 
AOM2 signals at ~ 80 MHz (c). These signals drive the acousto optic 
modulator AOMl and AOM2 and switch on and off the illumination 
and reference beams. 


gate duration St (that is much shorter than the vibration period 
Ta = 2'!TfujA) yields a broadening of the spectrum of about Ta/SI 
for m. In order to take quantitatively these effects into account, 
we have calculated the spectrum. The slowly varying complex field 
E{t) must be multiplied by a gate function H{t,td,St) of period 
Ta- Within the interval [OjTa] we have: 


H{t,td,St) = 0 

for 

0 < t < td — 5t/2 

= 1 

for 

td — 5t/2 <t <td + St/2 

= 0 

for 

td St/2 < t < Ta 


The slowly varying complex field signal E(t) becomes thus: 

E{t) = A (21) 

= A ^ A(m,0,5t,$) 

m 

Em = E 

where A{m,td,St) is the m**' Fourier component of the periodic 
function E{t)/E. We have calculated A{m,td,St,^) by Fourier 
transformation of E{t)/E, and we have plotted on Fig. [THl the 
Doppler spectrum \Em\‘^/\E\'^ = \A{m,td, St)]"^ for different gate 
times tc and for SI/Ta = 0.1. In Fig. [TH] (a), the gate coincides 
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Figure 18: Doppler spectrum \Em\‘^/\E\^ = \A{m,td,6t)\^ calculated for 
a vibration amplitude = 0.50, a gate width 6t = 0.1 Ta and a gate 
time tc = 0.25 (a), 0.35 (b) , 0.45 (c) and 0.50 Ta (d). 


with the maximum velocity {td = 0.5 Ta)- The Doppler spectrum 
is narrow and centered near m ~ — 4> = —50 rad. In figure [T5] 
(b) and (c) i.e. for gate time tc = 0.25 (a), 0.35 (b), the absolute 
velocity decreases, and the center of the Doppler spectrum moves 
towards m = 0. On the other hand, the Doppler spectrum broad¬ 
ens. In figure [TSl Id the gate coincides with the minimum velocity 
{td = 0.5 Ta)- The center of the spectrum is m = 0, and the 
broadening is maximum. 

To measure the instantaneous velocities of the object at time 
t = td, the gate is adjusted at td, and the holograms of all side¬ 
bands m are recorded. This can be done automatically by the 
computer, by adjusting the frequency u}lo to fulfil EqllHl The 
sideband images E.m{x,y) are then reconstructed for all m, and a 
3D cube of data (coordinate x, y and m) is obtained as illustrated 
by Fig. [TT] One can get velocity images by making cuts along 
X or y- Figure [TUI shows an examples of velocity images obtained 
at different times td of the vibration period. Here, the reed is 
oriented in the same direction than in Fig ll4l and the velocity im¬ 
ages of coordinates x, m are obtained with horizontal cuts in plane 
y = 256. 
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Figure 19: Successive velocities of the reed on a period Ta- These images 
are obtained by making a cut in the y = 256 plane of the 3D stack of 
reconstructed images, whose coordinates are x, y and m with m = -100 
to -1-100. Gate duration is 5t = 0.1 Ta- From (a) to (j), the gate time 
td is swept from = 0.25 Ta to = 1.15 Ta by 0.1 Ta step. The 
images are displayed in logarithmic scale for the sideband complex field 
intensity: \Em{x, y = 256)\‘^. 
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3.5 Sideband holography for small amplitude of 
vibration 

When the vibration amplitude <i> becomes small, the energy within 
sidebands decrease very rapidly with the sideband indexes m, and 
one has only to consider the carrier m = 0, and the two first 
sideband to = ±1. Time averaged holography [S] that detects 
the carrier field Eq is not efficient in detecting small amplitude 
vibration d), because Eq varies quadratically with d): 


Eo{<^) = E Jo{<^) 


( 22 ) 


On the other hand, sideband holography that is able to detect the 
two first sidebands fields E±i is much more sensitive, because E±i 
varies linearly with <i>: 


E±i{<^>) = E J±,{<^>) 

~ ±E $/2 


(23) 


This point has been noticed about 40 years ago by Ueda et al. m, 
who has made a comprehensive study of the signal-to-noise ratio 
(SNR) observed for classical (photographic film) sideband holog¬ 
raphy. The authors managed to observe vibration amplitudes of a 
few Angtroms, and linked the smallest detectable amplitude to the 
SNR in the absence of spurious effects. Later on, sub nanometric 
vibration amplitude measurements were achieved with sideband 
digital holography [miMiMiMiEi], and comparison with single 
point laser interferometry has been made |3H [30] . 

One must notice that the complex field amplitude E scattered 
by the sample without vibration depends strongly on the x y posi¬ 
tion. In a typical experiment, the sample rugosity is such that E is 
a speckle that is fully developed. The field E{x,y) is thus random 
in amplitude and phase from one speckle to the next. One cannot 
thus extract the phase of the vibration motion from a measure¬ 
ment made on single a sideband, for instance the sideband to = 1. 
To get the phase, one needs to measure the field components on 
several sidebands, for example the carrier field Eq and the first 
sideband Ei. For small amplitude 4), the carrier field ifo(4>) is 
very close to the field E = Eo{^ = 0) measured without vibration. 
By measuring both ifo(4>) and ifi(<i)) by sideband digital hologra¬ 
phy, one can then eliminate E in Eq. I^and Eq. [^getting by the 
way Jo(‘h)) •^i(‘&) and 4). This can be made by calculating either 
the ratio Ei/Eq [25] : 


El _ Ji($) 


(24) 


Eq Jo($) 


~ 4>/2 
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Figure 20: Sideband holography setup able to acquire the carrier and 
sideband signal simultaneously 


or the correlation Ei Eq\ 

\E, E*\ = \E\^ Ji(<i>) Jo(<i>) (25) 

~ \E\^^/2 

Both methods yield a signal (d)/2 or \E\‘^^/2) that is proportional 
to 

From a practical point of view, Eq and Ei are measured succes¬ 
sively in a standard sideband holography setup like the one of Fig 15] 
mnTj. A sequence of Umax frames with ujlo = loq — ujccdI^ 
is first recorded yielding Eq. A second sequence with ujlo = 
LlIq + CJA — ‘^ccdM is then recorded yielding Ei. 

It is also possible to record Eg and Ei simultaneously EH Hi 
EH EH- In that case, the RF signal that drives the second acousto 
optic modulator AOM2 is made of two frequency components (Fig. 
EH) at s-'^d 

^'aOM2 = ^aomi+locgd/^ (26) 

^AOM2 = ^AOMl -F WA — WcCn/4 

• The first component at ^aom 2 ^ whose weight is a, yields the 
carrier signal Eg, if the hologram El is calculated with demod¬ 
ulation on the -1 grating order i.e. with H = J2n=T 

• The second component at uiaom 2 ^ whose weight is /3, yields 

the sideband Ei with demodulation on the -Fl grating order 
i.e. with H = In- 
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Figure 21: Vibration amplitude Zmax (a) and phase ip (b) of a lamel- 
lophone, averaged over the first upper cantilever, versus excitation fre¬ 
quency a;A/(27r) (Hz). Insets: retrieved vibration amplitude and phase 
maps in the neighborhood of the resonance. The theoretical resonance 
line is in gray. 
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Figure 22: Amplitude maps of the out-of-plane vibration of a thin metal 
plate versus excitation frequency Holographic images at 40.1 

kHz (a), 61.7 kHz (b). Scalebar: 5 mm. 


Simultaneous detection of Eq and Ei, and calculation of d* by 
the ratio method of Eq. [M]has been used to study various vibrat¬ 
ing samples. Figure [2T] shows example of signals obtained with a 
lamellophone of a musical box [SU]. In that experiment, all the 
frequencies loa, ujccd, ^aomi, ^aom 2 and lo\om 2 are driven by 
numerical synthesizers. In this particular case, harmonics are ex¬ 
tracted by a temporal Fourier transform of a sequence of 8 raw 
interferograms. Here, optical sidebands of interest are downcon- 
verted to demodulated frequencies in the camera bandwidth. In 
order to avoid signal phase drifts from the measurement of one 
sequence to the next, AOM2 is driven using a phase ramp. As a 
consequence, the phases of Eq and Ei are well defined, and yield 
the phase (p of the mechanical motion. We have: 

p = Arg(£;i/Ao) + C (27) 

where C is a additive constant that depends of the relative phase 
of the synthesizers. Figure [2T] shows the averaged amplitude Zmax 
and phase p of the upper cantilever of the lamellophone, which is 
driven though its resonance frequency ~ 541.5 Hz. As expected, 
the phase p makes a jump of about tt when excitation frequency 
LUA crosses the resonance frequency. 

The method has been also applied to get full field images of 
surface acoustic waves as shown on Fig. mm- Comparison 
with single point Doppler vibrometry is reported in Fig. 1231 If the 
carrier and sideband components of the RF signal at 
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PZT voltage amplitude (V) 


Figure 23: Comparison of quantitative out-of-plane vibration ampli¬ 
tudes Zmax retrieved with the single-point laser vibrometer (□ symbols) 
and holographic vibrometer with /3 = a {o symbols) and P ~ 50a (-|- 
symbols). 
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Figure 24: (a,b) Reconstructed images of the correlation EiEq obtained 
with a cube of wood; (a) amplitude: \EiEq\ , and (b) phase: aigEiEQ. 
The correlation Ei(x,y)EQ(x,y) is averaged with a 2D gaussian blur 
filter of radius 4 pixels. The display is made in arbitrary linear scale, 
(c) 3D display of the phase: x,y,aTg EiEq 


uj\oM 2 have the same amplitude (i.e. if a = /3), some spurious 
detection of the carrier field Eq (which is much bigger) is observed 
when the detection is tuned on the grating order -1 (used to detect 
Eq). This spurious effect limits the detection sensitivity to Zmax ~ 
1 nm, as seen on Fig. [22] (o symbols). Much better results are 
obtained by reducing the weight of the carrier component (i.e. 
with P ':S> a). The detection sensitivity becomes Zmax 0.01 nm 
with P ~ 50q; as seen on Fig. (+ symbols). 

The ratio method of Eq. |23|used above is simple, since it gives 
directly 4). Nevertheless, the method is unstable for the points 
x,y oi the object where E is close to zero. Since 4) varies slowly 
with the location x, y, one can improve the accuracy by averaging 
4) over neighbor points. In that case, the ratio method is not 
optimal for noise, since all points x, y are weighted equally in the 
average, while the accuracy on the ratio depends strongly on the 
location x,y, since the accuracy is low when \E(x,y)\ is low. 

For this issue, the correlation method of Eqj^ is less simple, 
since it gives |£'p4), and not 4). It is nevertheless better for noise, 
since spatial averaging over neighbor points x,y is weighted by the 
scattered energy \E{x,y)p. 

Figure 1231 shows reconstructed images of the correlation EiEq 
obtained with a cube of wood (2 cm x 2 cm) vibrating at its res- 


30 



(a) 



Figure 25: (a) Clarinet reed with illumination beam focused in xo,yo- 

(b) Sideband m = 1 reconstructed image of the vibrating reed. The 
display is made in arbitrary log scale for the field intensity \Ei{x,y)\‘^. 


onance frequency uja = 21.43 kHz. Here, the fields Eq and Ei 
are measured sucessively. In order to get better SNR, the complex 
correlation signal Ei{x,y)EQ{x,y) is averaged over neighbor x,y 
points by using a 2D gaussian blur filter of radius 4 pixels. Fig¬ 
ure M shows the amplitude (a) and phase (b) of the the filtered 
correlation Ei{x,y)EQ{x,y). Figure [ 21 ] (c) displays the phase of 
the correlation in 3D. As seen, the opposite corners (upper left 
and bottom right for example) vibrate in phase, while the neigh¬ 
bor corners (upper left and upper right for example) vibrate in 
phase opposition. Note that the cube is excited in one of its cor¬ 
ner by a needle. This may explain why the opposite corners are 
not perfectly in phase in Fig. [211 (b,c). 

In order to increase the ability to analyze vibration of small 
amplitude, one can increase the illumination on the zone where 
the correlation averaging is performed. This is done by using a 
collimated illumination beam. Figure [21] shows an example. The 
illumination laser has been focused on point Xq , yo of the clarinet 
read to be studied (see Fig. |2^ (a) ). Most of the energy in 
the reconstructed image is thus located near Xo,yo (see Fig. | 25 ] 
(b) ). To measure the vibration amplitude at location xo,yo, we 
have calculated the averaged correlation (EiEq) and the averaged 
energy (|Aop): 



(28) 


Ji(4>)Jo(4>) {l/Nj,,,)Y,\E{x,y)\^ 
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Figure 26: Ratio {EiEq )/(vertical axis) as function of the reed 
excitation voltage Vpp in volt Units. The light grey curve is Ji(<h)/Jo(‘h) 

with d> = Zmaa:/(47rA). 


(iRoP) = {l/Np,.)J2\Eo{x,y)\^ (29) 

x,y 

= J2(<i>) (l/iVp„)^|R(x,y)|2 

x,y 


where y i® summation over the Npi^ pixels of the illumi¬ 
nated region located near xo,yo. In that region, the vibration 
amplitude is supposed to not depend on x and y. We get then: 


(iRoP) Jom 


(30) 


To verify the ability of measuring small vibration, the reed 
has been excited with lower and lower loudspeaker levels, the il¬ 
lumination laser being focused on the point xo,yo of Fig. [2^ (b), 
we want to study. A Sequence of Umax = 128 frames have been 
recorded for both carrier {Eq) and sideband (Ei). Four phase de¬ 
tection has been then made by calculating E[ with Unmax = 128 in 
Eq. [51 We have then calculated the reconstructed fields Eo{x,y) 
and Ei{x,y), the averaged correlation (EiEq) and averaged car¬ 
rier intensity (|ifoP) and the ratio (EiEq)/{\E o\'^). This ratio is 
plotted on Fig|55] {y axis) as a function of the peak to peak volt¬ 
age Vpp of the sinusoidal signal that excites the loudspeaker {x 
axis). The measured ratio varies linearly with Vpp. For low volt¬ 
age Vpp < 5 mV, the ratio reached a noise floor corresponding 
to {EiEq)/{\E o\'^) ~ 10“'*. Since this voltage Vpp is proportional 
to the vibration amplitude Zmax and thus to <i>, we have convert 
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the voltage Vpp in vibration amplitude Zmax- The conversion fac¬ 
tor {zmax = 0.02 nm for Vpp = 10 mV) is obtained by assuming 
that Ji(<i))/Jo(d>) fits the ratio {EiE^)/{\E q\^) that is measured, 
as predicted by Eq. [30] (see solid grey curve of FigI331l. We get 
here a noise floor noise corresponding to Zmax — 0.01 nm. It is 
about xlO lower than the limit A/3500 = 0.22 nm predicted by 
Ueda [26j . and xlO lower than in previous experiments at similar 
frequency I1Z1I2H] in the kHz range. Similar noise floor 0.01 nm 
has been obtained by the ratio method on Fig. but at higher 
vibration frequency a; ~ 40 kHz m- 

We have tried to make the experiment with more frames i.e. 
with Umax > 128, but this does not lower the noise floor. The 
remaining noise floor can be related to a spurious detection of the 
carrier held Eq when detection is tuned to detect Ei, or to some 
technical noise on the carrier signal that is not totally Alter off. 
Since the carrier and sideband fields Eq and Ei are within the 
same spatial mode, no space Altering can be applied to Alter off 
Eq. Here, the spurious carrier field Eq is Altered off in the time 
domain, while the local oscillator field E^o is Altered off both in 
time and frequency domain. 


4 Conclusion 

In this chapter we have presented the digital heterodyne hologra¬ 
phy technique that is able to fully control the amplitude, phase 
and frequency of both illumination and reference beams. Full au¬ 
tomatic data acquisition of the holographic signal can be made, 
and ultimate shot noise sensitivity can be reached. Heterodyne 
holography is an extremely versatile and powerful tool, in partic¬ 
ular when applied to vibration analysis. In that case, heterodyne 
holography is able to detect in wide field (i.e. in all points of the 
object 2D surface at the same time) the vibrating object signal at 
any optical sideband of rank m. Since the control of the inten¬ 
sity of the illumination and reference beams is fast, instantaneous 
measurements of the vibration signal, sensitive to the mechanical 
phase, can be made. For the measurement of large vibration am¬ 
plitudes, the possible ambiguity of measurements can be removed 
by making measurements at different sideband indexes m. For 
small vibration amplitudes, the mechanical phase can be obtained 
from measurements made on the carrier (m = 0) and on the first 
sideband (m =1). Moreover, extremely low vibration amplitudes, 
below 10 picometers, can be measured. 
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